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ABSTRACT 

We study the use of parallax microlensing to separate the effects of the mass function 
of dark massive halo objects (MHOs or 'machos') on the one hand and their spatial 
distribution and kinematics on the other. This disentanglement is supposed to allow a 
much better determination of the two than could be achieved entirely on the basis of the 
durations of events. We restrict our treatment to the same class of power-law spherical 
mode ls for the halo of MHOs studied in a previous paper (Markovic & Sommer-Larsen 
1997| ) . Whereas the duration-based error in the average MHO mass, p, = M/Mq exceeds 
(at N = 100 events) jj, by a factor of 2 or more, parallax microlensing remarkably 
brings it down to 15-20% of /2, regardless of the shape of the mass function. In addition, 
the slope a of the mass function, dn/dfi cx /i", can be inferred relatively accurately 
(<j a < 0.4) for a broader range, —3 < a < 0. The improvement in the inference of 
the halo structure is also significant: the index 7 of the density profile (p ~ R f ) 
can be obtained with the error a 1 < 0.4. While in a typical situation the errors for the 
parameters specifying the velocity dispersion profile are of about the same magnitude as 
the parameters, virtually all the uncertainty is 'concentrated' in linear combinations of 
the parameters that may have little influence on the profile and thus allow its reasonably 
accurate inference. 

Key words: Microlensing - Galactic halo - Macho mass function. 



1 INTRODUCTION AND OVERVIEW 

A statistical analysis by Alcock et al. (1996) of the 2-year 
microlensing data (6 or 8 events) obtained by the MA- 
CHO project in the direction of the Large Magellanic Cloud 
(LMC) indicated that the massive dark halo objects (MHOs 
or 'machos') responsible for the microlensing events could 
account for 30-100% of the total mass in the halo of our 
Galaxy. According to their analysis, typical (average) mass 
of the MHOs should lie in the range 0.1 - 0.6 M s . The more 
recent 4- year data (14 events, Axelrod 1997) yield similar 
ranges of the inferred quantities. Apart from the statisti- 
cal error due to the relatively small number of events, our 
ignorance regarding the structure of the halo of massive ob- 
jects (i.e., their spatial distribution and kinematics) leads to 
rather large uncertainties in the inferred masses. 

This last source of error is not likely to be extinguished 
if one relies only on the measurement of event durations 
T = -Re/^ii (Re is the Einstein radius, v n is the MHO's 



* Address after September 1st, 1997: Department of Physics, 
University of Illinois at Urbana-Champaign, 1110 W. Green St., 
Urbana, IL 61801, USA 



velocity orthogonal to the line of sight). Indeed, as shown by 
Markovic & Sommer-Larsen (1997, paper I), for the number 
of events N < 1000 the halo structure cannot be constrained 
sufficiently to allow a determination of the average mass 
p, = M/Mq of the MHOs to better than a factor of about 
2. Furthermore, paper I discussed only a limited class of 
spherical haloes; the results of Evans (1997), based on a far 
wider variety of halo models, imply that the range of JX (at 
virtually arbitrary N) could in principle extend from 0.1 to 
1. 

The duration T is, however, not the only relevant quan- 
tity that can be obtained from a microlensing event. For in- 
stance, photometric (Gould 1994a; Nemiroff & Wickramas- 
inghe 1994) or spectroscopic (Maoz & Gould 1994) methods 
have been proposed to measure the proper motion of the lens 
Vn/zD, where D is the Earth-source distance, z = D^/D 
and Dl is the Earth-lense distance. Another approach [dis- 
cussed by Grieger, Kayser & Refsdal (1986) in the context 
of quasar astronomy] is parallax microlensing, i.e., observing 
magnification through telescopes displaced from each other 
by about 1AU. More recently, Gould (1994b) studied and ad- 
vocated the use of parallax microlensing to obtain more in- 
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formation regarding the position and velocities of the lenses 
(and consequently reduce the uncertainty of their masses). 

The utility of the parallaxes stems primarily from the 
fact that the delay r between the maximal magnifications 
in the two detectors (one on the Earth and the other on a 
satellite in a heliocentric orbit) does not depend on the mass 
of the MHO crossing the two lines of sight to a source Addi- 
tional information is contained in the two maximal mag- 
nifications determined by the impact parameters Ui and 
U2 measured in units of the Einstein radius. Breaking a 
4-fold degeneracy (Gould 1994b; see also Section 2 of the 
present paper) by observing from a second satellite would 
allow us to obtain the so-called reduced transverse velocity 
v = v n /(l — z) of the MHO. A sufficiently large number of 
such measurements would then presumably put tight con- 
straints on the structure of the halo. However, even in ab- 
sence of a second satellite, the relative motion of the first 
satellite and the Earth could suffice to reduce the ambiguity 
to (at most) a 2-fold degeneracy in the direction of v for a 
majority of events (Gould 1995; Boutreux & Gould 1996). 
[According to Han & Gould (1995) this should (in the case 
of galactic bulge microlensing) permit measurement of in- 
dividual MHO masses to an accuracy of about 0.2 on the 
logarithmic scale.] 

In this paper we explore quantitatively the extent to 
which one could expect parallax microlensing to help con- 
strain the mass function of the MHOs as well as their spatial 
distribution and kinematics. The assumptions of the present 
paper are similar to those of paper I: for convenience we 
again limit ourselves to a class of spherical halo models (see 
section 3) described by a set of 5 parameters (the singular 
isothermal sphere is a particularly simple member of this 
class). On the other hand, the mass function is assumed to 
be a simple power law, dn/dfi oc /i Q characterised by three 
parameters (independent of the position in the halo): the 
average mass p,, slope a and range j3 on the logarithmic 
scale. 

Although lacking somewhat in generality, this frame- 
work will allow a straightforward application of the appara- 
tus of statistical parameter estimation: the errors of maxi- 
mum likelihood inference of the mass function and halo pa- 
rameters can be estimated from the sensitivity of the dis- 
tributions of directly measurable quantities to small shifts 
in the underlying parameters (see paper I and Section 5 of 
the p resent paper) . We call such estimates the Cramer er- 
rors (Cramer 1946j ). For simplicity we will study both the 
degenerate (with the full 4-fold degeneracy) parallax mi- 
crolensing, where the observable quantities are T, p = t/T, 
ui and lt2 and resolved (the 4-fold degeneracy completely 
removed) parallax microlensing with observables T, p and 
w — a 2 — p 2 , where a is the transversal distance between 
the two lines of sight in the lens plane measured in Einstein 
radii. 

In specific computations we at first adopt for the param- 
eters of the underlying halo model the values corresponding 
to the centrally condensed (7 = 3.4; p oc i? -7 ; R is the 
distance from the centre of the Galaxy) halo of blue hori- 
zontal branch field stars (BHBFS; see Section 3). Although 
one might speculate as to the relevance — or irrelevance — 
of this structure to the halo of MHOs (paper I), these val- 
ues are simply taken as a convenient starting point for our 



numerical experiment and the accuracy of their inference is 
estimated. In addition we briefly discuss the inference start- 
ing from the singular isothermal sphere (7 = 2; constant, 
isotropic velocity dispersion) as a model for the MHO halo. 

We find that parallax microlensing reduces the Cramer 
errors in jl from a factor of 2-10, characteristic of mea- 
surements of event durations, to only 15-20% (at TV = 100 
events) for a priori unknown halo structure parameters. This 
error is typical of inference under the (unrealistic) assump- 
tion that the halo model is accurately known a priori and is 
kept fixed in the maximum-likelihood fitting of the distribu- 
tion of measurable quantities. The improvement indeed re- 
sults from an effective disentanglement of the mass function 
from the halo structure. In addition, while parallax-based 
errors in parameters a and /3, specifying the shape of the 
mass function, are comparable to duration-based errors if a 
is sufficiently close to -1.5, the growth of the errors away 
from this value is strongly restrained by the parallaxes (see 
Fig. 10). 

Parallax microlensing also reduces by about two orders 
of magnitude the Cramer errors for the halo parameters. 
The power index 7 of the halo density profile is determined 
with the error <r 7 < 0.4 (at both 7 = 3.4 and 7 = 2), again 
from TV = 100 events. On the other hand, the errors in the 
parameters specifying the velocity dispersion profile are of 
roughly the same magnitude as the parameters themselves. 
This, however, does not necessarily mean that the profile 
will be poorly determined: as shown in Section 5, virtually 
all uncertainty (at 7 = 3.4) is due to the existence of a sin- 
gle linear combination of velocity dispersion parameters that 
is poorly constrained even by parallax microlensing. Indeed, 
small displacements in the parameter space along this vector 
cause a particularly weak change in the velocity profile. In 
addition, these changes tend to occur predominantly at large 
radii, where the microlensing rate is low. The peculiarities 
of the singular isothermal sphere, on the other hand, lead to 
the existence of three poorly constrained linear combinations 
of velocity dispersion parameters, none of them having a sig- 
nificant effect on the flat velocity profile. Remarkably, all the 
above conclusions are virtually independent of whether we 
use the degenerate or fully resolved parallaxes: the improve- 
ment due to resolving the degeneracy is modest. 

In Section 2 of this paper we derive expressions for dis- 
tribution functions of measurable quantities. These expres- 
sions are general and can be used with arbitrary halo mod- 
els. The specific class of halo models used in this paper is 
described in Section 3 (following a similar section in paper 
I). Section 4 deals with the morphology of the distribution 
functions derived in Section 2, while Section 5 explores their 
sensitivity to the underlying parameters and thus derives the 
Cramer errors of inference. Finally, section 6 contains basic 
conclusions of this paper along with a bit of speculation re- 
garding their more general validity. 



2 DIFFERENTIAL PARALLAX 
MICROLENSING RATES 

As observed from the Earth, a lens of mass M — /uM©, 
crossing the Earth-source line of sight at distance zD (0 < 
z < 1) from the Earth and with the impact parameter uiRe 
[Re = r Ey /JIy/z(l-z), r E = 2 v / GM Q D/c 2 = 3.2 x 10 9 km] 
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Figure 1. A parallax microlensing event projected onto the lens 
plane (orthogonal to the line of sight). All distances are given in 
units of the Einstein radius Re- 




will magnify the sta r by the maximum factor ^4. ma x = 
(u 2 + 2)/(u\y/uf + 4). On the other hand, a satellite will 
detect maximum magnification, determined by the impact 
parameter U2-Re relative to the satellite-source line of sight, 
with the time shift r from the moment of the Earth-observed 
maximum magnification. 

The geometric relations between the observable quan- 
tities, T = Re/vh (event duration as observed by both the 
Earth-based observer and the satellite), r, Ui and Ui, are 
easily derived from Fig. 1. If r is the component of the Earth- 
satellite vector orthogonal to the line of sight, its source- 
centered projection onto the plane (also orthogonal to the 
line of sight) of the lens is 



r(l 



Re 



r 1 



(1) 



(measured in units of the Einstein radius -Re). 

The points '1' and '2' in Fig. 1 denote intersections of 
the Earth-source and satellite-source lines of sight respec- 
tively with the lens plane. The lens' trajectory along the 
unit vector v = (cos <j>, sin <j>) (again, projected on the lens 
plane), shown as the solid arrow, crosses the 1-2 line at dis- 
tance uRe from the parallel line (dotted) drawn through the 
midpoint between 1 and 2. Consequently, 



r(l 



v. 



z) cos </> 1 
T 



T 

T 



(2) 



is an observable quantity. On the other hand, w = a sin <f> 
generally cannot be obtained unambiguously from u\ and 
U2 only; the 4-fold degeneracy is illustrated in Fig. ^. A 
secure way of breaking this degeneracy would be to use a 
second satellite [the line tangent to three circles of radii u\, 
U2 and M3 (the last measured from the second satellite) is 
unique]. However, as Gould (1995) has shown, the motion 
of the Earth and the single satellite relative to the line of 
sight should allow us to resolve the ambiguity — at least 
regarding the magnitude of w — in most cases. Nevertheless, 
in this paper we will discuss both the 'degenerate' (the 4- 
fold ambiguity unresolved) and the 'resolved' (w uniquely 
determined) parallax microlensing. 



Figure 2. Four possibilities corresponding to a pair of values 
tti > 0, U2 > in the case of degenerate parallax microlensing. 



The lens' rate of crossing near the line of sight per single 
source and a single (number density near the Sun n = 1) 
lens is 

d^-^— J DH(z)dz J RE{z)du J f n (v n , (f>)v^dv n d(j), 

(3) 

where H(z) [H(0) = 1] is the MHOs' halo density profile 
along the line of sight, f n [J f n {v n , 4>)v n dVn.d4> = 1] is the z- 
dependent, 2-dimensional distribution of velocities projected 
orthogonal to the line of sight and dn /dn [J (dn j d^£)dp, = 
1] is the z-independent mass function of the MHOs. 
Using 



(r/r E ) 



fia 2 + (r/rs) 2 ' 
rewriting the lens plane area element 

dpdw = adadd) = — f-r-rl dzdS, 
2 V da 2 ) 



(4) 



(5) 



and switching to integration over T (dv n /dT — —v B /T) we 
obtain 



x(T,p,w) 



dT 



dTdpdw 



= I d^DH(z)RE(z)f n (v n 



■tj 2(r/r E )V 
T [fia 2 + (r/r E ) 2 ] 



(6) 



where v n = TE^/Jly z(l — z) /T and tanci = w/p. 

By contrast with resolved parallaxes, where \{T,p, w) is 
of more immediate relevance, measuring u\ — \u + w/2\ and 
\u — w/2\ (see Fig. 1) is not sufficient to uniquely determine 
w in the case of degenerate parallaxes. Inserting 'dummy' 
integration J dui5(ui — \u + w/2\) J dii2<5(u2 — \u — w/2\) in 
the rate (0) and using the identity d(u+w/2) Ad(u — w/2) = 
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dw A du, we arrive at the expected result 



dTdpdu\du2 

= x(T,P, u i + wa) + x(T,p, -ui - u 2 ) 
+ x( T ,P, u i - u ?) + X( T ,P, - M i + "2), 



(7) 



expressing the differential rate in terms of variables accessi- 
ble to degenerate parallax microlensing detection. 

So far in this section we have ignored the question of the 
minimum amplification necessary for successful detection of 
a microlensing event. In particular, microlensing might pro- 
duce sufficient magnification only in the Earth- (or satellite-) 
based detector, i.e., the magnification in the other detec- 
tor (say '2') could be too small for a reliable determination 
of t and U2- To deal with this possibility, we will require 
Ui < «th, U2 < Uth for a detectable parallax (double, i.e. in 
both detectors) microlensing event, while single events will 
correspond to Mi < u th , 112 > u tn or iti > u th , U2 < it tn , 
where u t h is a certain threshold value. The detection rate of 
events (both single and double) in one (say '1') detector is 
thus 



P(T, Ul ) 



where 



fi (1) (T, Ul ) 



dTdm 

dU2 



du2 



ul single 

dTdui 



dp^{T,p,u 1 ,u 2 ) 
dpV(T,p,ui,u 2 ) +n (1) (T, Ul ), 



(8) 



dU2 



dp^(T,p, Ul ,u 2 ) 



(9) 



is the differential detection rate for single events, expressed 
in terms of the only available measurables, T and ui. Of 
course, fl w (T,u) = ft (2) (7» = Q.(T,u). The differential 
rate P(T,u), introduced in equation ([]), is 



P(T,u) = 2Dr% J d^^) 2 j\zHi 



z)[z(l- 



#/n("n,(/l). 



(10) 



If we assume (as we will in the present paper) that the 
MHO mass function can be well approximated by a simple 
power law 

dn 1 



du C (a) fx" 



_M 

Q + l ■ 



(11) 



where (3 = log 10 (/j ma x/Mmin), Mmax and /i min are the upper 
and lower bounds of the mass range, jj, = ^//^max^mm and 



Cp(a) = 



/31nl0 



a = — 1, 



)/2 



10 



-/3( a +l)/2] 



(12) 



the one-detector rate P(T, u) simplifies to 



P(T, u) = 2Dr% 



772(a+l) 

4 ± 



where y = n/T 2 and 



F(y) 



;/ I dz [z(l - 
>o 



wzHH 



'H(z) 



d<j) f n 



V a F(y)dy, 



(13) 



'■e 



\J z(l - z)y, 4> 



(14) 



In the rest of the paper we will use a 'composite' notion 
of microlensing event including double (parallax) events and 
single events detected only in detector '1' or '2'. The com- 
posite probability distribution function for degenerate par- 
allax microlensing can then be obtained by introducing the 
normalising constant A 

V(T,p,ui,u 2 ) = i*(T,p,^i,u 2 ), 



f7(T, u) 



Q(T, u) 



(15) 



so that 



dT 



dp 



du\ 



du 2 *(T,p, U!,U 2 ) 



+ 2 1 dT dun(T,u) = l. (16) 
Jo Jo 

In order to take account of the detection condition 
mi < Uth, U2 < uth in the case of resolved parallax mi- 
crolensing detection, we multiply x{T,p,w) by the range of 
u (see Fig. 1) for which the double event detection condition 
is satisfied 



X (T,p,w) (2ut h - \w\) x (T,p,w), 



(17) 



(\w\ < 2u t h), and thus obtain the differential detection rate 
of events characterised by the observables T, p and w. The 
relevant probability distribution for the detectable double 
events is then 



X(T, P ,W) = j(2u th - \w\) X (T,p,w), 



(18) 



where A has the same value as in the degenerate paral- 
lax case [one can show J™ th dui J Q th du 2 ^(T, p, Ui, u 2 ) = 
r 2 ^th ^ w f2u t h — \w\)x(T, p, w)], thus yielding the normali- 
sation 

dT I dp dwx(T,p,w) 



+ 2 



dT I duQ{T,u) = 1. 
'0 



(19) 



3 MODELS OF MHO DISTRIBUTION AND 
KINEMATICS 

In this paper we will consider a range of spherically sym- 
metric models of the massive halo objects' distribution 
and velocities. Probably the most commonly used model is 
the isothermal sphere with the velocity dispersion constant 
throughout the halo and the density profile which is well 
approximated by 
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Figure 3. Anisotropy parameter and velocity dispersion for 
the CS halo model as functions of the distance d from the Earth 



in the direction ot bMU 
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; cr r is given by tnc solid line, at = < 
ec paragraph fulluwnig ciinaLiuii (\24 



the dashed line. The straight solid lines correspond to the singular 
isothermal sphere, SIS (/3 = 0, cr = 156 km/s) . 



p{R) = p - 



■ + R 2 



(20) 



where a w 5 kpc is the 'core' radius and Rq = 8.5 kpc is 
the distance of the Sun from the galactic centre. Assum- 
ing that the total (luminous + dark matter) halo density 
is distributed according to expression (pcj), one obtains the 
observed (approximately) flat rotation curve for the galaxy. 

The MHO mass distribution, however, need not follow 
that of the total halo mass. We may, for instance, follow 
the hints provided by recent observations (Sommer-Larsen, 
Flynn & Christensen 1994, Sommer-Larsen et al. 1997) of 
the blue horizontal branch field stars (BHBFS) in the outer 
halo. These observations imply that the velocity dispersion 
changes from /3 = 1 — of / a i > (<7 t and a t are velocity dis- 
persions respectively in the radial and tangential direction 
relative to the Galactic centre) at smaller distances R from 
the centre of the Galaxy to /3 < at larger distances. The 
radial velocity dispersion is well described by the analytic 
fit 



1 1 
tan 

2 7T 



(21) 



where the best agreement with the observations is achieved 
with (j — 80kms _1 , <r+ = 145 km s -1 , r — 10.5 kpc and 
I = 5.5 kpc (these are the values used in paper I and the 
present one; more recent values, based on a larger sample of 
stars, are given in Sommer-Larsen et al. 1997). The BHBFS 
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Figure 4. Differential rate dT/dT (normalised) for the SIS (solid 
line) and CS (short-dashed) halo models with Sun's and LMC 
motion neglected. For the dotted line one assumes the CS model 
and takes into account the motion of the Sun and LMC. In all 
cases the mass function is a delta function centered on fi = 0.4. 



halo is close to spherical with the density that is well mod- 
eled by the power law p — p (Rq / R) 1 , where 7 ~ 3.4. 



The Jeans ' equation for spherical systems (Binney & 



Tremaine 1987) yields the tangential velocity dispersion 



2 1 T/ 2 



(i 



2 r da x 



(22) 



where V c = (— 7?d$/di?) 1//2 is the (roughly constant) rota- 
tion velocity. This tangential dispersion is smaller than in 
the case of an isothermal sphere (7 = 2, <T r = const.) with 
the same V c (see Fig. |). 

We will, following paper I, model the velocity distribu- 
tion by the Gaussian 



f(v r ,ve,v^) 



1 



1 



(2tt) 3 / 2 orcr; 



■ exp 



2 1 2 
v e +V4, 



(23) 

where <r r and a t are given by equations (|2l]) and ( |22[ ) for 
power-law density profiles. The relevant distribution [used, 
e.g., in the rate (H)] of velocities orthogonal to the line of 
sight is then 



fn(Vn,</)) 



■ exp 



1 / [V n COs((j) + (j)o) + Si] 

2 ~ 



+ ' 



[v n sin(0 + <j> ) 



(24) 

For the above equation we have introduced orthonormal 
vectors in the plane orthogonal to the line of sight: % is 
in the plane determined by the Sun, LMC and the Galac- 
tic centre (GC) and points in the general direction of GC; 
j = k x i, where k points along the Sun-LMC line of sight. 
Thus, Si and Sj are the corresponding components of the 
local (z-dependent) velocity of the line of sight relative to 
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Figure 5. Plots of x with the SIS halo (Sun's and LMC motion neglected) and a delta mass function (/i = 0.4) at T = 14.2 (top) and 
T = 84.7 (bottom) days. 



the galaxy 
a 



is the angle between the 



j — Ot and cr| 

~0 



: 5 of + sin 2 5 ( 



'2' axis and i, 
) cr," [sin 5 = (Rq / R) sin i 
and R 2 - = R 2 Q + (zDf ~2zDR Q cos t (t = 82° is the angle be- 
tween GC and LMC as observed from the Earth)] . Detailed 
derivations and numerical values are given in the Appendix. 



As in paper I, the halo model corresponding to BHBFS 
with the power-law density profile 7 = 3.4 and the disper- 
sion given by and Jeans' equation will be called the 
'concentrated sphere' ('CS'). More generally, we will assume 
for our study that the MHO halo can be described by a 
member of a class of models specified by five parameters: 7, 
(Jo, (7+ r and I. For instance, the model with 7 = 2 and con- 
stant velocity dispersion <r r = a t = K/\/2 = 156km/s (i.e., 



(Jo = V c /V2, <j+ 
sphere (SIS). 



0) is just the familiar singular isothermal 



4 DISTRIBUTION OF MEASURABLE 
QUANTITIES 

In order to understand basic features of the distribution 
functions x(T,p,w) [and consequently ^(T,p, Ui, u 2 )] and 
f2(T, u) we will at first limit ourselves to the relatively sim- 
ple case of the SIS halo, neglect the motion of the Sun (i.e., 
the detectors) and the LMC and assume that all MHOs have 
the same mass, [i = p, (see Fig. 4). 
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Figure 6. Plots of \ with the CS halo (Sun's and LMC motion taken into account) and a delta mass function (fj, = 0.4) at T = 15.1 
(top) and T = 81.4 (bottom) days. 



With the above assumptions 

. (l-z) 2 z 4 H(z) 
X (T,p, w) oc i >— A exp 



1 r|ju z(l - z) 



2 a 2 



rp2 



(25) 



where z(p,w) is given by equation (^). At sufficiently small 
a 2 = p 2 + w 2 , a MHO is near LMC, z w 1, 1-2 w (r E /r) 2 fia 2 
and we thus have 



X ~ a exp 



' 2 I ruT 



(26) 



Near the origin of the p, w plane \ grows at first as a 4 (from 
X — at a — 0) to a maximum value at oi = 2rcrT/r|/i 
(see the 'funnel' emerging from w — on the right of Fig. 5, 



which is a projection along the p-axis of the 3D plot on the 
left), which for the values of ths SIS model parameters and 
fx = 0.4, r = 2 AU takes on value oi = T/508 d. 

At large a, z w (r /t-e)' 2 / fJ-a 2 (close to the detector at 
2 = 0) and thus 



■ exp 



1 / r \ 2 1 
~2\~^t) a 2 



(27) 



This expression reaches a maximum at a-2 = (1/2 3 / 2 )r/aT 
7.9 d/T, and falls off as \ oc a -8 for a > a2. 
For T <T a , where 



1 r Ev /JI 
2V 1 a 



63 d, 



(28) 
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Figure 7. The fraction Q(T, u)/P(T, u) of events detected in one 
detector (i.e., with u < « t h = 1) that are not detected in the 
other detector. The upper panel gives the plots for the SIS halo 
model (without Sun's or LMC motion) and the bottom for the 
CS (with Sun's and LMC motion) model, both with a delta mass 
function, fi = 0.4. The lines are labeled with event durations T in 
days. The total (i.e., integrated over T) fraction of single events 
out of all events detected in a detector is 12% in the upper case 
and 20% in the bottom case: the CS model gives more events at 
low z, where i?g/r(l — z) is smaller. 



we have a\ < a,2 and an annular 'valley' will exist between 
the two maxima (see the top of Fig. 5) separating the inner 
'funnel' from the outer circular wall of the 'volcano.' The 
detected events would thus belong to two distinct classes 
distinguished by the magnitude of a. 

As we shift toward larger T's the 'valley' turns shallower 
and flattens out completely around T = T a . For T > T (see 




Figure 8. The lines tangential to both the circles of radius 
and uiig are drawn to find the sector of size 2(<j>i + fa) from 
which the events at u (< u t \^ = 1) in the right-hand side detector 
are not detectable by the one on the left-hand side. The distance 
between the centres of the circles is r(l — z). 



the bottom of Fig. 5) only one maximum at a = oi will 
persist. 

Introduction of the anisotropic velocity dispersion of the 
CS halo model will distort the circular contour lines of Fig. 5 
into elliptical ones. In addition, switching on the motion of 
the Sun and LMC (Earth's revolution around the Sun ne- 
glected; the Earth-satellite line lies in the ecliptic and chosen 
orthogonal to the line of sight; see the Appendix) will 'erode' 
the outer wall of the 'volcano' asymmetrically (Fig. 6): more 
MHOs will cross the 1-2 line of Fig. 1 in the upward direc- 
tion. The motion of the line of sight through the rest frame 
of the halo will also lead to a general shift of the differential 
event rate dT/dT toward shorter durations (Fig. 4). 

In Fig. 7 we plot for various durations T the fraction 
Q(T, u) /P(T, it) of events detected in one detector (say, u = 
U2 < Uth.) that are not detectable (iti > u t h) in the other 
one. At short T's (few days) the event rate is dominated by 
MHOs that cross the line of sight close to the Sun. For these 
MHOs Re = rEy^/Iy z(l — z) is smaller than the projected 
separation r(l — z) between the lines of sight with respect to 
the two detectors (a > 1). The elementary construction of 
Fig. 8 shows that of all MHOs passing at impact parameter 
uRe relative to detector 2, the fraction 



fl(T,u) 
P(T~u) 



H + 02 



— -(l-u)+cos — r(l + u) 

r(l — z) r(l — z) 



2 Re 



IT r(l — z) 37T 



Re 



r(l - z) 



(l + 3u 2 ) 



(29) 



will pass at impact parameters greater than Re with respect 
to detector 1 and thus not be detected if u t h = 1 (we here 
assume the 2-dimensional velocity distribution is isotropic as 
is the case for the SIS model without the Sun's and LMC's 
motion) . 

At R E < rfl - z) [z < 1/(1 + nrl/r 2 ) » (r/r E ) 2 /M = 
0.02] fraction (129) depends weakly on u\ the flatness of the 
corresponding curves can be observed in the top part of 
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o 1 

log[(r E /r)W(l-z))>/ 2 ] 

Figure 9. Normalised (i.e. of a constant area under the curve) 
dT/dTdz at T = 4.6 (solid), T = 11.2 (long dash), T = 166 
(dotted) and T = 361 days (short dash) for the SIS (without 
Sun's or LMC motion; top) and CS (with Sun's and LMC motion; 
bottom) halo models. In both cases the mass function is delta at 
H = 0.4. The 166d and 361d curves are virtually indistinguishable. 



Fig. 7. Naturally, the fi/P curves are depressed as T [and 
thus Rft/r(l — 2)] is increased. 

For MHOs passing at larger z > 0.02 (corresponding to 
typical times T « Re/V2ct > r/y/2a ss 15 d) the Einstein 
radius will be longer than the projected Earth-satellite dis- 
tance r(l — z). This means that an event of a small impact 
parameter u < 1 — r(l — z)/Re, with respect to one detec- 
tor, will inevitably be detected [Q(T, u)/P(T, u) = 0] by the 
other one. For u> 1 — r(l — z)/Re, the fraction detectable 
in only one detector rises approximately as 



fi(T, u) 



Re 



(30) 



1 -1 

r = — cos — ^ 

P{T,u) 7T r(l 

with increasing u and reaches 1/2 at u — Uth = 1. 

The correspondence between the event durations and 
typical z's is illustrated in Fig. 9. The differential rate 
dT/dTdz is relatively sharply peaked if plotted with respect 
to the coordinate (r E /r) [piz/(l — z)} 1 ? 2 — R E /r(l — z) — 1/a 
instead of z [the area element below a curve in Fig. 9 differs 
from (dT/dTdz)dz only by a constant factor]; we used this 
fact in the above discussion to relate a narrow range of z 
(or, rather, simply a single value of z) to each T. 



5 INFERENCE OF THE MHO MASS 

FUNCTION AND THE HALO STRUCTURE 

Even at first sight it seems obvious that the increased 
amount of information obtained from parallax microlens- 
ing measurements should allow a more reliable determina- 
tion of the structure of the halo and the MHO mass func- 
tion than just the measurements of event durations. Indeed, 
the measurable quantity pT — r(l — z) coscp/v n (equation 
^ involves only the kinematic properties of a MHO (i.e., 
excluding its mass). In the case of resolved parallaxes, an 
additional quantity, w, is measured and this gives us the 2- 
dimensional velocity projected on the observer's plane (the 
so-called reduced velocity; see. e.g., Gould 1994b) 

~1 r a r 1 . 

;(p,w). (31) 



r a 

T^I Vn = To 2 



Tp 2 + w 2 



Although the 4-fold ambiguity of the degenerate parallaxes 
reduces the quality of available information, both types of 
parallax microlensing would constrain the halo indepen- 
dently of the MHO masses; this should make it possible to 
separate the effects of the masses and determine the MHO 
mass function. 

In order to assess the information gain due to parallax 
microlensing we will investigate simulated maximum like- 
lihood inference of the five parameters (ci = 7, C2 = a , 
C3 = 0-+, C4 = r , C5 = I) of the halo model and the three 
parameters (c6 = u, cr = a, eg = P) of the mass function 
based on the normalised probability distributions x(T,p, w), 
Cl(T, u) (equation |l|) for resolved or *(T,p, Ui, Us), Cl(T, u) 
(equation hfl) for degenerate parallax measurables. 

Ideally, one would prefer to obtain an answer through 
a series of Monte-Carlo simulations, but the enormous nu- 
merical task of computing the above 3-dimensional distri- 
bution x f° r many points in the 8-dimensional parameter 
space compels us to seek more economical alternatives. For 
a sufficiently large number iV of detected events, the average 
\ 1 / 2 

■j, C^i '• ) ] 2 \ of the parameters (cji is the 

'real' value of a parameter) can be shown (see, e.g., paper 
I) to approach asymptotically the Cramer limit 



(32) 



where C M „ is the inverse of the information matrix 



N 
+ 2N 



dTdPduidu 2 



dTdu 



(33) 



(, M denotes derivative with respect to the fi'th parameter) 
for degenerate parallaxes and by analogy for the resolved 
parallaxes. 

In our numerical experiment we assume that the pa- 
rameters of the concentrated sphere (CS) describe accu- 
rately the MHO halo. The Cramer errors of the inference 
of the mass-function parameters for iV = 100 are shown in 
Fig. 10 as functions of the mass-function slope a. The re- 
solved parallax errors obtained with the assumption that 
both the halo structure and the mass function parameters 
were unknown (and thus variable in the maximum likelihood 
fitting) are given as solid lines. For comparison, we also show 
the Cramer-limit errors (dotted lines) for measurements of 
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Figure 10. The Cramer-limit errors with an unknown (variable) 
halo from resolved parallaxes (solid lines) and event durations 
only (dotted). For a fixed CS halo model, the errors are shown 
as dashed lines (resolved parallaxes) and dot-dashed (event du- 
rations only) lines. All errors are given for p, = 0.4, (3 = 2 and 
N = 100 events. 



event durations only. In addition, we plot the resolved par- 
allax errors (dashed lines) and the errors based on event du- 
rations only (dot-dash) assuming that the halo parameters 
are known precisely (and thus not varied for the maximum 
likelihood fit). In all computations we take into account the 
motion of the Sun and LMC, assume r = 2 and that the 
Earth-satellite segment is in the ecliptic and orthogonal to 
the Earth-LMC line of sight (see the Appendix). 

The most striking feature of these plots is a significant 
reduction of the error in p, brought about by parallax mea- 
surements. As shown in paper I, by changing the parameters 



of the halo model one can match closely event duration dis- 
tribution curves corresponding to widely different average 
masses. This is reflected in the large errors in p, if the halo 
parameters are allowed to vary in addition to the parame- 
ters of the mass function. The extra information provided 
by parallaxes effectively allows us to constrain the halo as to 
bring down to the values comparable to those of inference 
with a fixed (i.e., 'known') halo model. (Notice that the gain 
in accuracy due to parallaxes for a known halo structure is 
modest by comparison.) 

In paper I we have concluded that for a certain range of 
a (very roughly, — 2 <, a <, 0) the inference of a and /3 is rel- 
atively weakly affected by the uncertainty of the halo struc- 
ture. This again is manifested in errors of the same (small) 
order near a = —1.5 in all four cases shown in Fig. 10. How- 
ever, away from a = —1.5, a a and op grow large for event 
duration-based inference. Here again, the parallaxes restrain 
this growth significantly (especially for a a ) and keep the er- 
rors down in the range characteristic of the fixed-halo in- 
ference. Interestingly, degenerate parallax errors a a and 
ct/3 are only a few percent smaller than the resolved parallax 
values (and we omit the corresponding plots). 

As suggested above, the significant improvement in the 
accuracy of the mass function determination should be as- 
cribed to the disentanglement of the halo structure from 
the MHO masses. This disentanglement is obvious from the 
plots in Fig. 11 of the correlation matrix 



C in/ 



(34) 



The off-diagonal spikes 'coupling' the halo with the mass 
function are noticeably less pronounced in the parallax case 
(right). [Notice, e.g., the strong correlation between p, and 
7 in the duration based inference (left) and its significant 
drop in the parallax case (right).] 

Further comparison of the two plots of Fig. 11 also 
shows a significant reduction in the correlation among some 
of the halo parameters due to the parallax measurements. 
Does this imply that the halo structure itself could be in- 
ferable? Figure 12 indeed displays a remarkable suppression 
of the Cramer errors of the halo parameters' inference: the 
(resolved) parallax errors at least are of the same order of 
magnitude as the parameters themselves or smaller In par- 
ticular, the density profile index 7 can be determined rather 
accurately. [The difference between the degenerate and re- 
solved parallax errors is more pronounced for the halo model 
parameters than for the mass function; still, the degenerate 
parallax errors are larger than the resolved parallax ones by 
at most 20%. The measurement of p and the information 
contained in ui, 112 and Q(T,u) are thus sufficient to con- 
strain the halo even if the determination of w is subject to 
the 4-fold ambiguity.] 



t As discussed in paper I, the Cramer errors can serve as reliable 
estimates of actual errors only if they are small in comparison with 
the corresponding parameters that one is trying to infer. Typi- 
cally, we would expect the nonlinear dependence of the distribu- 
tion function of the observables on the underlying parameters to 
make the actual errors smaller that the Cramer estimates when 
these are relatively large. In this case the actual errors would fall 
slower than TV -1 / 2 and approach the Cramer limit from below. 
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Figure 11. Correlation matrices for Cramer errors of inference based on event durations only (left) and resolved parallaxes (right). The 
underlying halo model is CS and the parameters of the mass function function are p = 0.4, a = —1.5 and (3 = 2. 



O-(l) 


= 0.033, 


V(i) 


CT (2) 


= 0.795, 


V(2) 


f(3) 


= 0.126, 


V (3) 


CT(4) 


= 0.383, 


V(4) 


O-(B) 


= 3.316, 


V(5) 


f(6) 


= 0.182, 


V(6) 


0-(7) 


= 0.102, 


V(7) 


C(8) 


= 0.081, 


V (8) 



( 0.7667, 0.2005, 

( 0.0362, 0.7289, 

( -0.4744, 0.2243, 

( 0.0511, 0.2676, 

( -0.0158, -0.5050, 

( 0.3327, -0.1966, 

( -0.2675, 0.1085, 

( -0.0270, 0.0343, 



0.5791, 


0.1562, 


-0.0838, 


-0.5037, 


0.4946, 


0.1677, 


-0.3075, 


0.8129, 


0.3601, 


-0.0903, 


-0.3563, 


-0.1221, 


0.2431, 


0.0859, 


0.0510, 


0.0499, 



-0.1037, 


0.0233, 


0.4541, 


0.0229, 


-0.0685, 


0.3439, 


0.4111, 


0.0083, 


0.7789, 


-0.0143, 


0.0442, 


0.8117, 


-0.0271, 


0.4702, 


0.0048, 


0.0186, 



0.0149, 


-0.0258 ) 


-0.0048, 


0.0035 ) 


-0.5733, 


0.0175 ) 


-0.0103, 


-0.0338 ) 


-0.0019, 


-0.0001 ) 


-0.2121, 


0.0540 ) 


0.7796, 


-0.1434 ) 


0.1351, 


0.9871 ) 



Table 1. Eigcn-errors (for N = 100 events) and the corresponding eigcn- vectors of the relative error matrix for the CS halo and fj, = 0.4, 
a = -1.5 and (5 = 2. 



The residual but stil considerable correlations among 
the errors of the halo parameters indicate that some com- 
binations (functions) of them (together with the mass func- 
tion parameters) may be possible to determine with a sig- 
nificantly higher accuracy. In order to investigate this possi- 
bility we examine the eigen-values and eigen- vectors of the 
relative error matrix 

^iMltHte- 1 ))' (351 

i.e., those normalised vectors V( p ) {V( p )^(p) = 1j using Ein- 
stein's convention for summation over the repeated index fj.) 
that satisfy the eigen-equation 

A^Vft^Vfo. (36) 

These eigenvectors correspond to the linear expansions of 
the desired functions of the parameters in the vicinity of 
c^°' ; they indicate the mutually independent (uncorrelated) 
infinitesimal displacements in the parameter space caused 
by statistical errors of parallax microlensing detections. In 
a certain sense, they are the eigen-modes of the halo model 
(plus the mass function) as viewed through the 'instrument' 
of parallax microlensing. 



We illustrate the above point by a specific example. 
Table 1 contains the square roots ('eigen-errors') of the 
eigen-values along with the eigen- vectors of the relative er- 
ror matrix, assuming again the CS halo model and fl = 0.4, 
a = —1.5 and (3 = 2 (same as for Fig. 11). We immedi- 
ately notice that the eigen-errors for the eigen-'modes' arc 
small {(Jip) < 1) except for a single mode (a (5) = 3.316). 
This, the 5'th mode 'mixes' rather strongly C2 = a , c 5 = I 
and to a lesser degree C3 = a + (see the corresponding com- 
ponents of V(5)). Its mass function components, V^, VZ^ 
and V(|), are very small which explains why its large asso- 
ciated eigen-error does not lead to large errors of the mass- 
function parameters. On the other hand, although there is 
an eigen- vector, V( 3 ), that strongly mixes the halo and the 
mass mass function parameters, its relatively small eigen er- 
ror, <T(3) = 0.126 contributes little to the errors of the mass 
function parameters and thus gives at most moderate cor- 
relations between these and the halo-model errors. 

In Fig. 13 we compare changes in the halo model caused 
by a displacement in the direction of the vector Vru of 
a small eigen-error with a displacement along Vr^y Notice 
that a much larger displacement (e = 0.5) is needed along 
V( 5 ) to cause a discernible effect. In addition, these changes 
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Figure 12. The Cramer-limit errors (at p, = 0.4, (i = 2 and 
N = 100) of the halo-model parameters from resolved parallaxes 
(solid lines) and event durations only (dotted lines). 



tend to be concentrated at larger distances, where the event 
rate (recall 7 = 3.4) is smaller; microlensing is, not surpris- 
ingly, relatively insensitive to such displacements along V( 5 ) 
in the parameter space. As a consequence, the large associ- 
ated eigen-error, (7( B ), is the leading culprit for the relatively 
large errors of C2 = a , c 5 = / and C3 = o + evident in Fig. 12. 
By contrast, parallax microlensing constrains fairly well the 
displacements along the other directions. In particular, the 




20 30 40 
d [kpc] 



Figure 13. Comparison of the lowest-error halo mode, Vni, at 
e = 0.05 (left) and the highest-error, V( 5 ) at e = 0.5 (right) halo 
mode. In the bottom plots <r r of the CS model is given as long 
dashed lines while its 'shifts' corresponding to ^(1) or eV(5j 
are given as dot-long dashed lines. The ctj = at lines are short 
dashed (CS) and dot-short dashed (shifted from CS). The ex; lines 
are solid (CS) and dotted (shifted). In all cases p, = 0.4, a = —1.5 
and /3 = 2. 



eigen-modes closest to the parameters of the mass function 
(V(6), V(7) and V(s)) all have rather small eigen-errors. 

The coefficient e = 0.5 used in Fig. 13 for the displace- 
ment along V(5) is considerably smaller than the associated 
eigen-error 075) = 3.3 (see Table 1); the relatively small e 
allows us to stay near the 'linear' range of shifts in the struc- 
ture of the halo (the range of the validity of the Cramer er- 
rors). At e = 3.3 the structure would be changed radically, 
far beyond what one might expect on the basis of small 
linear displacements. This likely implies that actual errors 
would be smaller than the Cramer-limit estimate. The re- 
sults of such strongly non-linear shifts can be studied only 
by means of Monte-carlo simulations which are beyond the 
scope of this paper. 

The CS halo discussed so far can be regarded as a typ- 
ical member of our class of spherical models. On the other 
hand, SIS is peculiar in that the parameters r and / can be 
arbitrary as long as 0+ = 0. Might this 'degeneracy' lead to 
a large uncertainty in the inference of the halo structure? 

To explore this issue we assume 7 = 2, a = 156 km 
0+ = 10 km, r = 10.5 kpc and I = 5.5 kpc, only slightly 
deviating from SIS. The root-quadratic dependence of a T 
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CT (3) V (3) 


°"(4) V (4) 


0-(5)V {5) 


7 


2.02 


1.81 


1.77 


a [km/s] 


170 


116 


15 


a+ [km/s] 


81 


111 


272 


r [kpc] 


-1.7 x 10 3 


-1.1 x 10 3 


-2.5 x 10 3 


I [kpc] 


-0.6 x 10 3 


-0.5 x 10 3 


4.3 x 10 3 


cr r [km/s] 


186 


122 


158 


at [km/s] 


155 


160 


164 



Table 2. Halo structure parameters and radial and tangential 
velocity dispersions after displacements cr^j V«) in the parameter 
space. 



on a and o + (equation Ell) suggests we use the squares 
c 2 = al = 2.4 x 10 4 km7s 2 and c 3 = a\ = 100 km 2 /s 2 
among the parameters whose inference we estimate in the 
Cramer limit. 

At iV = 100 events (we take p, = 0.4, a = — 2 and 
(3 = 2) the resolved parallax errors for the density profile, 
o\ = 0.36, and the mass-function parameters, oe = 0.08, 
(77 = 0.20, as = 0.27, are again small, while the errors of 
the velocity-dispersion parameters are much larger: 02 = 
2.9 x 10 4 km 2 /s 2 , a 3 = 8.1 x 10 4 km 2 /s 2 , 04 = 3.4 x 10 3 kpc 
and <T5 = 4.7 x 10 3 kpc (the degenerate parallax errors are, 
again, only insignificantly larger). 

The diagonalisation procedure used above for the CS 
case gives for the eigenvectors mixing the velocity-dispersion 
parameters <T( 2 ) = 0.30, cr(3) = 207, <J( 4 ) = 183 and <7( 5 ) = 
1100. Although the three eigen-errors 07 3 ), au\ and arg) 
seem inordinately large, the corresponding shifts crnVm in 
the parameter space do not produce dramatic changes in 
the structure of the halo. Table 2 gives the values of the 
halo parameters after these shifts, as well as the radial and 
tangential velocity dispersions. Not surprisingly, in all three 
cases the very large magnitudes of r a and I ensure position- 
independent velocity dispersion profiles and through very 
different combinations of a and cr+ lead to a r and cr t rela- 
tively close to those of SIS. This 'degenerate' (almost) SIS 
structure shows even more clearly than the above CS-based 
example that large uncertainties in the values of the halo 
parameters are not incompatible with a good grasp of the 
structure of the halo: parallax microlensing can probe the 
halo while ignoring the vagaries of our parametrisation of 
its structure. 



Moreover, some combinations of parameters describing the 
halo structure can be inferred rather accurately. Virtually all 
the uncertainty regarding the halo model is then localised 
in a few (precisely one in the CS-based example discussed 
in Section 5) 'eigen-modes' of the halo, i.e., those displace- 
ments in the halo parameter space that are left 'loose' by 
parallax microlensing while allowing other displacements to 
be independently (and tightly) constrained. These poorly 
constrained modes correspond to particularly small changes 
in the actual structure of the halo as given by density and 
velocity dispersion profiles. Unfortunately, their large asso- 
ciated eigen-errors may push us into a non-linear regime 
(inadequately charted by the Cramer limit) of deviations in 
the halo structure that can be properly explored only by 
time-consuming Monte-Carlo simulations. 

Although obtained on the basis of a limited class of 
halo models, the above conclusions should be relevant in a 
broader context. The halo structure/mass function disentan- 
glement (and consequently accurate mass determination) as 
well as the inference of some properties of the halo struc- 
ture itself should result from the significant enhancement of 
information (as elaborated in Section 5) due to parallax mi- 
crolensing even if one allows for a much wider range of halo 
structures. For instance, in the more general case of non- 
spherical haloes, one may hope to constrain the density and 
velocity dispersions along the line of sight well enough to 
determine the mass function. The 'uncertain' modes would 
then describe possibly large changes in the halo structure 
away from the line of sight. Only parallax microlensing ob- 
servations in several directions would presumably suffice to 
infer the overall structure of the halo. These conjectures need 
to be tested on other, more realistic or better dynamically 
founded models (e.g., those of Evans 1994) of the halo than 
the ones discussed in the present paper. 
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Appendix: Position and motion of the Earth 
relative to the Galaxy and the Large Magellanic 
Cloud 

This appendix contains formulae (with derivations) giving 
positions and proper motions in the rest frame of the Galaxy 
of points along the Earth-Large Magellanic Cloud line of 
sight. 

We denote by x, y and z unit vectors directed, respec- 
tively, toward the Galactic centre (GC), tangentially along 
the rotation of the 'local standard of rest' (LSR) around the 
galactic centre, and toward the north Galactic pole. If the 
velocity of a source relative to the Sun is given in terms of 
Galactic-coordinate components w ra( j, Vi and vi, the conver- 
sion to the x, y, z components is 

Vj = (cos b cos I u ra d — sin b cos Ivb — sin I vi ) x 
+ (cos b sin I v ra ,d — sin b sin I Vb + cos lvi)y 
+ (sin b w ra( j + cos bvb) z (Al) 

in the rest frame of the Sun. From the 1990 data of B. Jones 
(as quoted in Greist 1991) the proper motion of the LMC is 



Wrad 
Vb 
VI 



250 ± 5km/s 
335 ±62 km/s 
-31 ±62km/s. 



Disregarding the error bars 
b = -32°.8, I = 281°, 
55kpc (0.160, -0.825, -0.542); 
vUc = (44, -390, 146) km/s. 



(A2) 

and using D — 55kpc, 
[thus r L MC = -Dk = 

k ■ k = 1] we arrive at 
Bearing in mind that the 



Sun's velocity relative to the Galaxy is vq = (0, 220, 0) + 
(9,11,16) = (9, 231, 16) km/s, where the first and second 
terms are the velocity of LSR relative to the Galaxy and 
that of the Sun relative to LSR, the velocity of the LMC 
relative to the Galaxy is 



vlmc = v + vlmc = (53, -159, 162) km/s. 



(A3) 



[After completing the calculations described in this paper, 
the author learned about more recent values of Jones, Kle- 
mola & Lin (1994) who give v LM c = (60 ± 59,-155 ± 
25, 144 ± 51) km/s. On the other hand, Kroupa & Bastian 
(1997) give an estimate based, among other sources, on the 
Hipparcos data: v LM c = (41L44, -200L31, 169L37) km/s. 
These values do not differ sufficiently from velocity ( Ac ) to 
lead to a significant change in the results reported in this 
paper.] 



Along the Sun-LMC line of sight we can introduce 
auxiliary coordinates defined by the following orthonormal 
system: as above, k is directed along the line of sight, i 
lies in the Sun-GC-LMC plane and poins in the direction 
of the GC and j = k X %. Given the data listed above, 
i = (0.987,0.134,0.088) and j = (0,-0.549,0.836). Each 
point on the line of sight moves relative to the galaxy with 
velocity s(z) = Vq(1 — z) +vlmcz whose components along 
I and j are 



V Qi (l — Z) + DLMC.Z 
V Qj (1 — z) + VtMCjZ, 



(A4) 



where vq { = 41 km/s, VhMCi = 45 km/s, VQ j — —113 km/s 
and Vi,MOj = 223 km/s. 

In the above we have neglected the motion (revolu- 
tion) of the Earth (detector '1') and the satellite (detector 
'2') around the Sun: one should in principle add the term 
u> x r® (or the corresponding term for the satellite) to vq. 
To avoid (at this stage) unnecessary complication we will 
keep disregarding the revolution around the Sun. For spe- 
cific computations in this paper we assume that the 'l'-'2' 
line points along the unit vector f directed along ill x k. 
Since in galactic coordinates ub — (—0.095,0.862,0.498), 
(derived from data in Mihalas & Binney 1981) we find 
f = (-0.652, 0.324, -0.685) and thus r • i = -0.660 = cos 4> a 
and f • j = —0.750 = sin</> D corresponding to (f> = 229°. 



